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oo Abstract 

We investigate the holographic fermions over a gravitational lattice background with a rather 
CN 
t^J- low temperature. Since the rotation symmetry is broken on the plane, the lattice effects change 

o 

C*") the shape of the Fermi surface within the first Brillouin zone from a circle to an ellipse. When the 



Fermi surface intersects with the Brillouin zone boundary, the band structure with a band gap is 



^ observed through a numerical analysis. We construct a lattice model sourced by a scalar field as 

well as an ionic lattice model without the scalar field. In both cases we find the similar physical 
results. 
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I. INTRODUCTION 

As a powerful tool of understanding a variety of strongly coupled condensed matter sys- 
tems, AdS/CMT has been entering an era from building purely theoretical models with 
highly symmetric configuration to modelling more realistic condensed matter system with 
less symmetry. In particular, instead of working within the probe limit, very recently the au- 
thors in jl] and [2] have constructed some spatially inhomogeneous but periodic gravitational 
backgrounds by fully solving the coupled partial differential equations numerically with the 
Einstein-DeTurck method. Such gravitational backgrounds, by holography, correspond to 
the boundary systems in the presence of a lattice, a key ingredient in the condensed matter 
systems. Thus with such gravitational backgrounds, one can explore various lattice effects 
in the holographic investigation of condensed matter systems. Actually for the transport 
coefficients such as optical conductivity and thermoelectric conductivity, it has been shown 
that not only does the presence of a lattice result in the Drude peak at low frequencies 
but also induces a new intermediate scaling regime in which a robust power law behavior is 
found with respect to the frequency [H [2]. Remarkably such a result is in striking agreement 
with the experiments on the cuprates, including the superconducting phase [3]. 

The purpose of this paper is to investigate how such a holographic lattice affects the 
Fermi surface by putting the Dirac field in such bulk gravitational backgrounds. To achieve 
this, in the next section we first build up the holographic framework to extract the spectral 
function for fermions living on the boundary by solving the bulk Dirac equation in periodic 
backgrounds. Then we shall present a numerical construction of two kinds of ultra cold 
holographic lattices using the Einstein-DeTurck method in Section III. After that, the nu- 
merical results for the Fermi surface in the presence of such holographic lattices are detailed 
in Section IV. We shall conclude with some discussions as well as further directions in the 
final section. 

II. HOLOGRAPHIC SETUP FOR THE BULK DIRAC FIELD IN PERIODIC 
BACKGROUNDS 

Start from the bulk action for a Dirac field with mass m and charge q 

S D = i [ d 4 xV^( (T a V a - m) C, (1) 



in the static background as general as 

ds 2 = —g tt (x, z)dt 2 + g zz (x, z)dz 2 + g xx (x, z)dx 2 + g yy (x, z)dy 2 + 2g xz (x, z)dxdz, 
A = A t (x,z)dt. (2) 

Here T a = (e M ) a r M with (e fl ) a a set of orthogonal normal vector bases and T M Gamma 
matrices. In addition, T> a is the covariant derivative, i.e., 

V a = d a + ^(u lxl/ ) a T^-iqA a , (3) 

where Y^ u = ^[T^,T U ], and (w M ^) a are the spin connection 1-forms given by 

(uv)a = (e M ) 6 V a (e,) 6 . (4) 

With the following orthogonal normal vector bases 
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Thus the Dirac equation 

T a V a ( - m = (7) 



can be written as 



— ^— r r 3 d z ( + -^=r°(d t - i q A t )( + -Lr 1 + 9x ' 

9x X gzz-g xz y/gtt wfe v gxx{g X xgzz - g xz 



-r 3 d x ( 



H-=r 2 9 y c + 



1 

'9yy 



9 X gtt OxQyy 



+ 



^gtt^[g*~x ^g yy ^g~x~x A(g xx g zz - g 2 xz ) 

1 



gxxO x g zz g X z\gxzO x g xx ^gxxQxgx 



3/2 



kgJx (g X xg Z z - g. 



2 ^ 

XZ) 



W9xx(gxxg zz - glz 

gxz 



F) Of) J_ i,xz f) _l_ 9xxOz9tt gxzO x gtt gxxO z gyy gxzO x gyy \ -p3 

o z g xx ^o x g xz -\- o x g xx -\- + 

gxx gtt gyy 



(8) 



To proceed, let us first make a transformation ( = (gttgxxgyy) A J~ . Then the above equation 
turns out to be 
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Next expanding J 7 as J 7 = -F(x, z)e 8U, *+ zfc > x * ; one can have 



A 3 r 3 F + A r°F - Air x F - A 2 r 2 F + mF = 0, 



(10) 
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Now if we choose our gamma matrices as 
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and split the 4-component spinor into two 2-component spinors as F = (Fi,^) , then the 
Dirac equation becomes 

A 3 * | - ma 3 (8) I * ] + A a 2 ® f * ± i/\ x a l ®\ \ - i^a 1 ®l 2 = 0. (13) 
^2/ \F 2 J \F 2 ) \F 2 J \F t J 

Note that this is a coupled equation between F\ and F 2 as opposed to the case in [6] , where 
iq and F 2 are decoupled. Furthermore, by the decomposition 



with a = 1,2, the Dirac equation (13) can be expressed as 
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Suppose our background field is periodic along the x direction with the periodicity c. 
Then by the Bloch theorem, the solution to our Dirac equation can always be expanded as 
follows 

A*0£j2)\ _ >r^ / A a ,n{Z) J inKx QgN 

B a (x,z)l n=0,±l,±2,.» \Ba,n{z) 

with K = ^. 

c 

Below we will restrict ourselves into the following periodic background, i.e. 
ds 2 = ±-{-(l- z )P( z )Q a ( x , z )dt 2 + Qzz{x > z)dz2 



z 2 " y ' y ^" v ' ' P{z){l-z) 

+Q xx (x, z) [dx + z 2 Q xz (x, z)dz] 2 + Q yy (x, z)dy 2 }, 

A = (1- z)ip(x,z)dt (19) 



Here 

P{z) = 1 + Z + Z i_dl_ (20) 

with the other variables regular from the horizon z = 1 all the way to the conformal boundary 
z = 0. In particular, we require at the horizon 

Q tt (x,l) = Q zz (x,l) (21) 

and at the conformal boundary 

Qtt(x, 0) = Qzz(x, 0) = Q xx (x, 0) = Q yy (x, 0) = 1, Q xz (x, 0) = 0, ij,{x, 0) = fi(x). (22) 

Whence such a geometry asymptotes AdS with the curvature radius L — 1 and the temper- 
ature given by 

pm = 6-j4 (23) 

4vr 8tt V ; 

With the above background ansatz, our Dirac equation gives rise to 

at the horizon. In order to obtain the retarded Green function on the boundary by holog- 
raphy the independent ingoing boundary condition should be imposed at the horizon, i.e., 

A \ / 1 

1 - z) ^t (25) 

B / \ — i ' 

for each a and n with the others turned off. On the other hand, near the AdS boundary, 
our Dirac equation reduces to 

(zd z - ma 3 ) ® [ Fl ' n 1=0. (26) 

Whence the solution can be asymptotically expanded near the AdS boundary as 

F a , n ^a a , n z m ( X J + b a , n z- m I ° | . (27) 



Holography tells us that the retarded Green function can be obtained by the following 
relation 

0-a,n((3,l) — G a>n . a i >n rb a ' !n '((3,l), (28) 



where a a%n (P, 1) and b a ^ n (f3 } I) are the asymptotic expansion coefficients in (27) of the solution 
to the Dirac equation evolving from the ingoing boundary condition with the only (/3, 1) mode 
turned on. 

In what follows, for simplicity but without loss of generality, we shall work solely with 
the case of m = 0. 

III. NUMERICAL CONSTRUCTION OF ULTRA COLD HOLOGRAPHIC LAT- 
TICES 

A holographic lattice background can be constructed at least in two ways. One way is 
to introduce a neutral scalar field with the periodic boundary conditions along the spatial 
direction. We would like to refer to this kind of holographic lattice as the scalar lattice. 
The other kind of holographic lattice is induced directly by a periodic chemical potential 
to the gauge field on the boundary, which is referred to as the ionic lattice. In the latter 
case, one does not need a scalar field. Both of these two kinds of holographic lattice can 
be numerically constructed by the Einstein-Deturck method, which has been detailed in 
[U HI E]. For our purpose, namely in order to see the lattice effect on the Fermi surface, we 
will follow such a routine to construct these two kinds of holographic lattice backgrounds 
with a rather cold temperature as follows, respectively. 

A. The scalar lattice 

To construct an ultra cold scalar lattice, let us consider the following gravitational action 
including Maxwell field and a neutral scalar field, i.e., 



S = ir n I (Pxy/=g 



R + — - -F ab F ab - 2V a $V a $ - 2M 2 $ 2 
L 2 2 



(29) 



where L is the AdS radius as before, and M is the mass of the scalar field. In what follows, 
we shall set L = 1 and M 2 = —2. Then the equations of motion can be derived from the 



above action as 



Rab + 3gab ~ 2 (V a $V fe $ - <$> 2 g ab ) - (F ac F b c - ^F cd F cd ) = 0, (30) 

v a n = o, (31) 

□$ + 2$ = 0. (32) 



Now with the ansatz in Eq.(19) for our holographic lattice, the asymptotical behavior of 
scalar field goes like 

$ = Z( j) = z[<\> x + z<f)2 + 0(z 2 )} (33) 

near the AdS boundary. Furthermore, the periodic structure can actually be induced by 
setting the source 

4>i(x) = A cos(k x), (34) 

and keeping the chemical potential fixed as 

H(x) = /i. (35) 

For our purpose we would like to construct an ultra cold lattice whose temperature is 



controlled by the parameter Hi in Eq.(23). This can be achieved by solving the coupled 
Einstein-Maxwell-scalar equations numerically using the Einstein-DeTurck method. Our 
numerical method is changing the partial differential equations into non-linear algebraic 
equations by the standard pseudospectral collocation approximation, and then solving them 
by employing a Newton-Raphson method. As an example, we show a solution to these 
equations in Figjljwith A = 1.5, k = 2, jjl = /ii = 2.35, which corresponds to a lattice with 
T/fj, = 0.0081. Note that different from the scalar field whose period is given by 27r//c , the 
corresponding solutions of all the components of the metric as well as the gauge field have a 
period of 7t/ko along the x direction, which comes essentially from the fact that in our first 
equation of motion the stress tensor is quadratic in the scalar field $. Since the scalar field 
does not appear in the Dirac equation, we claim that the lattice constant that the probe 



Dirac field feel is 7r/k . Therefore, the parameter K in Eq.(18) is given by K = 2k . 

It is worthwhile to point out that during the course of numerical analysis the lower the 
temperature is, the harder monitoring the accuracy is. In order to guarantee the convergence 
of our method at the given temperature, we are required to demonstrate the decaying ten- 
dency of the charge discrepancy An which is defined as A^ = |1 — Qn/Qn+i\ with Qn the 
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FIG. 1: We show Q zz , Q xz , ip and </> for k = 2, A = 1.5, \i = 2.35 and T/fi = 0.0081. 



charge on the boundary when the number of grid points takes N. As shown in in Figj2j such 
a decay is exponential, implying our results are exponential convergence with the increasing 
of the grid points. Moreover, we have also checked the behavior of the DeTurck vector field 
£ a which is defined in pQ and found that for our solution £ a £ a is smaller than 10~ 10 , ensuring 
that our numerics is leading to an Einstein solution rather than a Ricci soliton. 
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FIG. 2: On the left we show the boundary charge density p as a function of x, which can be read 
off by expanding -0 = H + (/■* ~~ z 9 )- 2 + @( z2 )- On the right we show the charge discrepancy An as a 
function of the number of grid points N, where the boundary charge is defined as Q = L ° dxp. 
The vertical scale is logarithmic, and the data is well fit by an exponential decay: log(Ajy) = 
-2.96- 0.20 N. 



B. The ionic lattice 



Compared to the scalar lattice, the ionic lattice can be generated simply by turning off the 



scalar field in (29) and imposing the spatially varying boundary condition for the chemical 
potential as 

ip(x, 0) = /i[l + Aocos(kox)]. (36) 

As an example, we illustrate a solution of Q xz and ip in Fig|3] with ko = 2, Aq = 0.1, 
/i = fii = 2.3, which corresponds to T/fi = 0.01. We stress that this type of solutions are 
different from those for the scalar lattice. Namely the chemical potential is periodic rather 
than a constant in the scalar lattice. As a result, the boundary charge density in the ionic 
lattice is expected to vary more dramatically than that in the scalar lattice. This can easily 
be seen in our plot for the charge density in Fig|3j In this sense, one will also expect that 
the effect onto the Fermi surface due to the ionic lattice should be much stronger than that 
due to the scalar lattice. We shall show in the subsequent section that this is actually the 
case. 
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FIG. 3: We show Q xz and ip for ko = 2, Aq = 0.1, /i = 2.3 and X '/ ' \x = 0.01 in the top and 
the corresponding boundary charge density as well as the tendency of charge discrepancy in the 
bottom. 

IV. NUMERICAL RESULTS FOR THE FERMI SURFACE 

The existence of Fermi surfaces for holographic fermionic liquids has been testified in 
various circumstances. We refer to [7] and [8] for a recent review. In this paper we focus 
on reporting the following two new results relevant to the Fermi surface when a lattice 
is introduced for the gravitational background. One is the shape of Fermi surface. In 
general one expects that the shape of Fermi surface is not circular any more since the 
rotation symmetry is broken in the presence of the lattice along the x direction. But how 
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the Fermi surface would change appears obscure since the Dirac equation becomes very 
complicated in this case. Even in the case in which the lattice effect can be treated by 
conventional condensed matter approach, there is no some kind of universal result on the 
shape of Fermi surface. Instead the shape will depend on the specific behavior of periodic 
potential. Nevertheless to our surprise, our numerical analysis gives us a very simple and 
elegant answer to this issue. That is, when our holographic Fermi surface is located within 
the first Brillouin zone, the shape is always an ellipse irrespective of the specific value of 
parameters in question. The other is the emergence of band gap at the intersection of the 
Fermi surface with the Brillouin zone boundary. This phenomenon coincides perfectly with 
the familiar lattice effect as one expects. Now let us demonstrate our numerical results in 
detail. 

To proceed, we would like to remark on how to identify the Fermi surface in our current 
setting. Note that we are working with holographic fermionic liquid in the presence of 
lattice at very low but non-zero temperature, where the Fermi surface is somewhat ill defined, 
because the concept of Fermi surface can only be defined at zero temperature with translation 
symmetry unbroken. However, as pointed out in |9j, the ARPES experiment is actually blind 
to the lattice by smearing it into a continuum somehow and the measured spectral function 
can be captured by the imaginary part of diagonal components of retarded Green function, 
namely A(u,k x — k\ + nK,k y = k 2 ) = Im(G 1]n . ln + G 2i n;2n)- Furthermore, taking into 
account that our ultra cold lattice only smears the Fermi surface in a negligible way, we 
can locate the position of Fermi surface by searching the peak of A(u, k x , k y ) with the tiny 
frequency u in the momentum space. Such an identification of Fermi surface is also described 
in [10] and similar to the operational definition given in [11]. 

With this in mind, we show in Figjl] an example with the 3d plot of A(oo, k) for the 
fixed k y = 0. A sharp peak occurs near u = 0, indicating a Fermi surface with the Fermi 
momentum kp = 2.634. Similarly, we may locate the positions of the Fermi surface for other 
values of k y such that the shape of the Fermi surface can be plotted in the momentum space. 
We illustrate our results in FigJ5]for the case of ko = 2, A = 2, fx = 2.35, T//J, = 0.0081 and 
q = 1.3, which corresponds to a Fermi surface located within the first Brillouin zone, namely 
kp < K/2 = 2. Although the shape of Fermi surface appears like a circle in plotting, our 
data clearly tells us that it is not a circle any more. As a matter of fact, it can be precisely 
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FIG. 4: We show the 3d plot of A(u,k x ,0) for k = 2, A = 1.5, (i = 2.35, T//J, = 0.0081 and 
q = 1.7. A sharp peak near u = implies a Fermi surface with the Fermi momentum kp = 2.634. 



fit by an equation of ellipse as follows 



k 2 k 2 
a 2 b 2 



y -1 



(37) 



with a = 1.8991 and b = 1.8511. We find this fitting has a very high accuracy, which can 
be seen from our error bar analysis presented on the right side of Figj5j Remarkably such a 
result is universal in the sense that our elliptical Fermi surface is robust against the values 
of parameters in question except that the longer axis as well as the shorter axis is varied 
as it should be. Now let us see how our elliptical Fermi surface varies with our relevant 
parameters. As such, we introduce two quantities, namely, the difference between the longer 
axis and shorter axis d = a — b as well as the eccentricity e = \/a 2 — b 2 /a. The relevant 
results are listed in Table [Tj, [TTJ and III Obviously, we observe the following behaviors 



As we increase the amplitude of our periodic source A Q , both d and e are increased. 
As we lower the temperature by increasing /ii, both d and e are increased. 
As we increase the charge q, the Fermi surface is enlarged with e suppressed. 



We remark that above phenomena can be observed in the ionic lattice background as 
well. Interestingly, we find that in ionic lattice case, the longer axis of the ellipse changes 
from x to y direction. 
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FIG. 5: We show the shape of Fermi surface on the left and the error bar fit by the ellipse in the 
polar coordinate on the right for ko = 2, Aq = 2, fj, = 2.35, T/fi = 0.0081 and q = 1.3, where the 
Fermi momenta kp < 2, namely the Fermi surface is located within the first Brillouin zone. 



A Q 


0.5 
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1.5 


2 


2.5 


d 


0.004 


0.013 


0.029 


0.048 


0.069 


e 


0.06414 


0.1187 


0.1742 


0.2242 


0.2675 



TABLE I: The variation of d and e with the amplitude Aq, where we have fixed the other parameters 
as k = 2, fj, = 2.35, T/fi = 0.0081 and q = 1.3. 

When the charge q is large enough, the Fermi surface will go beyond the first Brillouin 
zone with kp > K/2, Now let us turn to such a situation by first demonstrating the 
band structure of Fermi surface in Fig|6] for the scalar lattice. Note that the Fermi surface 
exhibits a periodic structure along the k x direction as it should be the case guaranteed by 
Bloch theorem. On the other hand, it appears that our Fermi surface does not show the 
band gap structure at the intersection of the Fermi surface with the Brillouin zone boundary, 
which is at odds with the familiar lattice effect. In order to see if this is really the case, 
we zoom in the Fermi surface precisely at the Brillouin zone boundary. The corresponding 
result is plotted in FigjT} To our pleasure, the band gap shows up and becomes large with 
the increase of the amplitude of periodic source. So to see a clearer band gap structure, it is 
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Ml 
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2.1 


2.2 


2.3 
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0.019 


0.021 


0.024 


0.027 


0.031 


e 


0.1487 


0.1546 


0.1610 


0.1696 


0.1792 



TABLE II: The variation of d and e with the temperature which is controlled by [Xi, where we 
have fixed the other parameters as ko = 2, A$ = 1.5, \x = 2.35 and q = 1.3. 



q 


0.5 


0.7 


0.9 


1.1 


1.3 


d 


0.013 


0.018 


0.022 


0.026 


0.029 


e 


0.2479 


0.2258 


0.2010 


0.1871 


0.1742 



TABLE III: The variation of d and e with the charge q, where we have fixed the other parameters 
as ko = 2, Aq = 1.5 and fi = fi\ = 2.35. 

better to move on to the ionic lattice in which the lattice effects should be much stronger as 
we discussed in previous section. As shown in Figj8]for the ionic lattice, the gap is evidently 
observed as one expects. From this figure we also notice that the band structure becomes 
richer when the Fermi surface intersects with more than two Brillouin zones. 




FIG. 6: Band structure for q = 1.7,/cq = 2, Aq = 1.5 and // = (jl\ = 2.35(kp > 2). 
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FIG. 7: There are two peaks for the imaginary part of the Green function at k x = 2(BZ boundary). 
The left one is for Aq = 1.5 while the right one for Aq = 2, with q = 1.7, ho = 2, and \i = \i\ = 2.35 
fixed. 
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FIG. 8: Band structure in ionic lattice model for kp > K/2{q = 1.7, Aq = 0.3, /U = fX\ = 2.3) and 
Uf > K(q = 1.7, Aq = 0.1, /x = \i\ = 2.3). Note that for the left one we have set ko = 4 such that 
iC = 4, while for the right one ko = 2 such that K = 2. 
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V. CONCLUSION 

We have investigated the lattice effect on the Fermi surface by putting the bulk Dirac field 
in the ultra cold holographic lattices, where both the background equations and the probe 
Dirac field are solved by pseudo-spectral method numerically. One interesting result is that 
the Fermi surface is always modified by the lattice from a circle to an ellipse. Suppose that 
such a Fermi surface is formed by some kind of exotic free fermionic quasi-particles, then the 
elliptical shape means that the effective mass along the x direction in scalar lattice (or the 
y direction in ionic lattice) becomes more massive by some kind of renormalization effects 
due to the presence of the lattice. Anyway, such a universal behavior of holographic Fermi 
surface begs for a deeper understanding. In addition, we holographically reproduce the band 
gap structure at the intersection of the Fermi surface and the Brillouin zone boundary, which 
is kind of well known lattice effect on the Fermi surface in condensed matter theory. 

Note that such a band gap structure is also obtained in |9J, we would like to take such 
an opportunity to compare our work with that made in [H|. The lattice considered in that 
paper is sort of weak ionic lattice, where the back reaction of the periodic chemical potential 
to the bulk geometry is ignored. The advantage of such a weak potential limit is two fold. 
One is that such a lattice can be put at zero temperature, and the other is that the analytic 
technique developed in [6] can be borrowed heavily to the relevant perturbation calculation. 
The disadvantage is also obvious. First, the amplitude of spatially varying chemical potential 
must be small enough, otherwise the perturbation calculation will break down. Second, the 
neglect of the back reaction may lead to the loss of some significant physics. For example, the 
elliptical Fermi surface we find does not appear in the perturbation calculation 1 . Third, when 
such a perturbation calculation is made to higher order or for the non-diagonal components 
of retarded Green function, the analysis will become much involved. Compared with this, 
the full retarded Green function can be obtained by our numerics once and for all. In this 
sense, our paper and jH] are complementary to each other. 

We conclude with various issues worthy of further investigation. First, numerically harder 
though it is, it is apparently better to construct the zero temperature holographic lattices 
for our probe Dirac field to propagate. Second, pertaining to the STM experiment, it is 



1 Of course, such a Fermi surface may show up when one goes to higher order in the perturbation calculation. 
After all, the perturbation calculation is made only to second order in [9]. 
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important for us to extract the retarded Green function in the position space, which can 
be achieved in the following two ways. One is to Fourier transform our resultant Green 
function in the momentum space to the position space. The other is to work directly in the 
position space by imposing Dirac delta source on the AdS boundary. We hope to address 
these issues in the near future. 
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